Coexistence of superconductivity and ferromagnetism in ferromagnetic metals 
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In this paper we address the question of coexistence of superconductivity and ferromagnetism in 
the high temperature superconductor RuSr2GdCu208-i. Using a field theoretical approach we 
study a one-fermion effective model of a ferromagnetic superconductor in which the quasiparticles 
responsible for the ferromagnetism form the Cooper pairs as well. We discuss the physical features 
which are different in this model and the standard BCS model and consider their experimental 
consequences. 



Recently an itinerant ferromagnet undergoing a high 
temperature superconducting transition was discovered 
in the copper oxide compound RuSr2GdCu208-5 (Ru- 
1212) and the experimental studies have revealed 

that the ferromagnetic state exists even below the su- 
perconducting transition. This prompts the interesting 
question of the possible many-body itinerant fermionic 
systems supporting both types of broken symmetry j^] . 

The search for ferromagnetic superconductors goes 
back to the sixties when superconducting materials with 
magnetic impurities were studied Q. The research in 
this direction has led to the works of Larkin and Ovchin- 
nikov 0, and Fulde and Ferrell js) (LOFF), who studied 
a simple model of effective field theory of superconduct- 
ing fermions coupled to magnetic impurities and they 
described the phase diagram of such a system. 

The Ru-1212 has similarities with the above materials, 
but is generically different because the fermions in the 
Ru02 layers (with mostly ferromagnetic fermions) are 
hybridized with the fermions in the Cu02 layers (with 
mostly superconducting fermions). Motivated by the 
above experimental observations we consider an effective 
two band Hamiltonian with ferromagnetic exchange in- 
teraction in one of the bands. One can perform a unitary 
transformation to diagonalize the kinetic term, including 
the tunneling between the two bands. As a result of 
this hybridization one obtains two new fermionic fields. 
Both fermions have ferromagnetic and pairing interac- 
tions (either strong/weak or weak/strong) and there are 
cross terms. The approximation of this paper involves 
neglecting the cross terms. Then one has two decoupled 
one-fermion problems. The one of interest in this paper 
is for the mostly superconducting fermions, because their 
exchange coupling is much smaller than the exchange 
coupling in the Ru layer, but could be comparable to 
Tc. 

In this paper for the first time the self consistent equa- 
tions for the superconducting gap and the magnetization 
are solved simultaneously in the mean field limit. We 
study the one fermion model of a ferromagnetic super- 
conductor leaving the two fermion model with the cross 



terms for further investigation. This case is relevant to 
the doped Ru02 layers in Ru-1212, where the magnetism 
becomes itinerant and the Ru02 layers participate in the 
transport properties of the material. 
Our model Hamiltonian is 

H-fiN^ J d'rciir) (^^V^ - c.(f) (1) 
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where Ca{r) are the spin a fermion fields, S ~ ^c\,TacriCai 
is the spin field, are the Pauli matrices, and /i is the 
chemical potential. The exchange interaction is ferro- 
magnetic (J > 0) and the four fermion interaction is at- 
tractive (g > 0). This is the simplest model which leads 
to the coexistence of ferromagnetism and superconduc- 
tivity. 

The partition function of the model can be written as a 
functional integral over the Grassmann fields c(t, r) and 
c{T,r) [0. We introduce a real vector field a1(t, r) us- 
ing the Hubbard-Stratonovich transformation of the ex- 
change term and a complex scalar field /(t, r) using a 
Hubbard-Stratonovich transformation of the second term 
in Eq.(l). The vector field describes the fiuctuations of 
the magnetization, while the complex scalar field the su- 
perconducting fluctuations. Performing the Gaussian in- 
tegral over the fermionic fields we obtain the partition 
function of the model as an integral over M, / and / 
which we calculate using the steepest descent around the 
mean field solutions M — (0,0, M) and A = g < f >. 
Here M = — < > defines the magnetization of the 
system. The mean field equations are 

where F^f f is the Free energy of a theory with the effec- 
tive Hamiltonian 

^^ff H [44t^pT + 44i^Pl + ^c_p-xCp-| + Act^cLp-^ 
p 



1 



2m* 



JM 



'P 2m* 



JM 



(3) 



Here the fermionic effective Hamiltonian H^f f is obtained 
after the Hubbard- Stratonovich transformations and set- 
ting the fields at their mean field values. 

Next we diagonalize the effective Hamiltonian using a 
Bogoliubov transformation. After the transformation the 
new dispersion relations are 
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where tp = — fi. Then the mean field equations take 
the form 
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where and are the momentum distribution func- 
tion of the Bogoliubov fermions and we have assumed 
that |A| ^ 0. 

From Eqn.(|^) one sees that for A/ > (the conven- 
tion that we will use here) E'^ > for all momenta p 
and therefore for T = 0, rip — 0. For E^ there are two 
possibilities. When JM < 2|A|, E^ < for all p and 
therefore — 1 for all p. Substitution of this in Eqn.(5) 
leads to M = 0. Therefore the only solution of the mean 
field equations which allows for the coexistence of fer- 
romagnetism and superconductivity is in the case when 
JM > 2|A| which we will assume. Then the equation 
E^ — has two solutions: 
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The dispersion of the /? fermion is positive when pp < 
p < and is negative in the complementary interval. 
The p dependence of E^ is depicted in Fig.l. With this 
in mind the Eq.(5) and (6) at T = have the form 
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It is difficult to solve analytically these equations, how- 
ever when JM is greater, but close to 2 A, p J is approx- 
imately equal to pp and therefore M is small as follows 
from Eqn.(8). In this case one can expand the r.h.s. of 



Eqn. (7) in the small parameter ^y(JM)^^^^A\A\^ obtain- 
ing 
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where pf — ^J2^jLm* . Substitution of these expressions 
in Eqn. (8) shows that in this approximation the magne- 
tization is linear in |A|, namely 



lAI 



(11) 



where r = Jm*pF / ^t^^ and this expression is valid for 
large r (i.e. MJ - 2|A| 0+). 

As in the standard BCS theory of superconductivity, 
the pairing of the quasiparticles occurs in the vicinity of 
p F , which must include the interval between p^ and . 
Then the integration in the first integral on the r.h.s. of 
Eqn. (9) is limited to a shell of width 2A, i.e. 
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Here + A > and — A < pp. 




FIG. 1. E!^ as a function of p. 

Substitution of the approximate expressions for 
from Eqn. (10) in the second term in the r.h.s. of Eqn. (12) 
and using the expression for the magnetization from 
Eqn. (11) we see that this term is independent of |A| and 
only leads to a renormalization of the cutoff A. In that 
approximation the solution is 
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When the magnetization increases the domain of inte- 
gration in the second integral on the r.h.s. of Eqn.(12) 
can exceed the size of the domain around pp in which the 
pairing occurs and which is the integration domain in the 
first integral of the same equation. In that case the sec- 
ond integral dominates and this leads to the absence of 
solutions with a finite gap. Taking the limiting case when 
the two integration domains are equal, i.e. pp -I- A = 
and pp — A = where are the values of the momenta 
from Eqn.(lO) with A = 0, we obtain the critical value 
of the magnetization 
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above which the superconductivity disappears besides 
the existence of an attractive four fermion interaction. 

Next we calculate the distribution functions and nj; 
of the spin up and spin down quasiparticles. In terms 
of the distribution functions of the Bogoliubov fermions 
these momentum distribution functions are 



(16) 



where and are the coefficients in the Bogoliubov 
transformation. They are independent of the magnetiza- 
tion and have the same form as in the BCS theory. 

At zero temperature is zero and — 0{pp — p) + 
0{p~p+). Then the spin up and spin down quasiparticles 
have the following momentum distribution functions 
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The functions are depicted on Fig. 2. 

The appearance of the Fermi surfaces of the Bogoli- 
ubov fermion (3 is unexpected in the superconducting 
phase, but it is a necessary condition for the existence of 
itinerant ferromagnetism. Therefore in the case of coex- 
istence of superconductivity and ferromagnetism caused 
by the same quasiparticles the existence of the two Fermi 
surfaces is a generic property of this state. These Fermi 
surfaces are reflected in the spin up and spin down mo- 
mentum distribution functions as well as in the anoma- 
lous Green's functions. It is easy to show that the anoma- 
lous Green's function, 

T{t - t',p) = - (Tcx(r, -p)ct(T',p)) , (19) 

in the case r = t' is 



T{Q,p) = 



(20) 



when < p < pp and p > p'^ and is zero when the 
momentum p is between the two Fermi surfaces p~^ and 
P > Pp- 
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FIG. 2. The zero temperature momentum distribution 
functions for spin up and spin down fermion. 

The existence of the Fermi surfaces, leads to different 
thermodynamic properties of the system, compared to 
the standard BCS theory. The specific heat has a linear 
temperature dependence at low temperatures as opposed 
to the exponential decrease of the specific heat in the BCS 
theory 
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Here 
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is the sum of the density of states on the two Fermi sur- 
faces of the Bogoliubov fermion f3. When the magneti- 
zation is small, from Eqns.(lO) and (11) follow that the 
density of states increases with r as N{0) -^rpp, 
as opposed to the case of ordinary weak ferromagnets, 
where the density of states is A^(0) = i^pp in this limit. 
Hence, the specific heat is large even at very low temper- 
atures. In the case of a superconductor in an external 
magnetic field |^ although there are gapless fermionic 
excitations the specific heat is not linear as opposed to 
our case. This can also be contrasted with some of the 
unconventional superconductors which have power law 
dependence of the specific heat on the temperature, de- 
pending on the nodal structure of the gap function. 
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Another consequence of the existence of the Fermi sur- 
faces is the existence of paramagnons which describes the 
longitudinal spin fluctuations They exist in ferro- 

magnetic normal metals and in our theory they survive 
even in the ferromagnetic superconducting phase. Their 
propagator is given by 



where a, b, and S are constants. The constant 
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defines the analytical properties of the paramagnon and 
is different from zero because of the existence of the Fermi 
surfaces. The constant 5 is 



5 = l-iN{0) 



(25) 



and 6 is a positive constant. As we mentioned earlier, 
the density of states, Eqn.(22), increases as the magne- 
tization, M decreases and therefore, at small, but finite 
value of the magnetization, M = Mq, S becomes zero, 
as opposed to the weak ferromagnetic metals where S 
becomes zero at zero magnetization. In the case of co- 
existence of the superconductivity and ferromagnetism 
the superconductivity prevents the magnetization from 
becoming arbitrarilly small, because when the magneti- 
zation is smaller than the critical value Mq, 6 is negative 
and the paramagnon fluctuations lead to an instability of 
that phase. The superconducting phase, with zero mag- 
netization (BCS like regime) the spin fluctuations of the 
paramagnon type are absent. 

Recently, a band structure calculation was performed 
by Pickett et al. Q and they have studied the origin of 
the superconducting state in the Ru-1212 compound. In 
our paper we considered the possibility of the coexistence 
of ferromagnetism and superconductivity and the physi- 
cal features of such a system. We arrived at a system of 
self consistent equations for the magnetization and the 
superconducting gap, and solved analytically these equa- 
tions at small magnetizations. The solutions with coexis- 
tence of superconductivity and ferromagnetism describe 
Bogoliubov fermions one of which has two Fermi surfaces. 
Therefore the spin up and spin down quasiparticles have 
two Fermi surfaces each. The thermodynamic properties 
of the coexistence phase are different from the standard 
BCS theory. The specific heat has a linear temperature 
dependence as in normal ferromagnetic metals, but in- 
creases anomalously at small magnetizations. These re- 
sults are obtained in a mean field approximation, but 
they are generic for the coexistence state and can be used 
as a starting point for a calculations beyond mean field. 
In our model the quantum critical point is dressed, i.e. 
the superconducting state occurs at zero magnetization, 
because the superconducting gap is generated not by the 



spin fluctuations, but by some other means. This is to 
be contrasted with the theory of spin fluctuations me- 
diated pairing in weak ferromagnetic metals [|2| where 
the quantum critical point is naked and the supercon- 
ducting ferromagnetic critical temperatures go to zero at 
the quantum critical point. In this paper we have con- 
sidered only uniform states. However, periodic solutions 
(like the LOFF state in the magnetic impurity case) will 
likely exist for certain regions of parameters, but possibly 
involving periodic magnetic structures as well as modu- 
lated superconducting order parameter. 
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